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ABSTRACT 

The Tremaine- Weinberg equations are solved for a pattern speed that is allowed to vary with radius. 
The solution method transforms an integral equation for the pattern speed to a least squares problem 
with well established procedures for statistical analysis. The method applied to the Hi spiral pattern 
of the barred, grand-design galaxy NGC 1365 produced convincing evidence for a radial dependence 
in the pattern speed. The pattern speed behaves approximately as 1/r, and is very similar to the 
material speed. There are no clear indications of corotation or Lindblad resonances. Tests show that 
the results are not selection biased, and that the method is not measuring the material speed. Other 
methods of solving the Tremaine- Weinberg equations for shearing patterns were found to produce 
results in agreement with those obtained using the current method. Previous estimates that relied on 
the assumptions of the density- wave interpretation of spiral structure are inconsistent with the results 
obtained using the current method. The results are consistent with spiral structure theories that allow 
for shearing patterns, and contradict fundamental assumptions in the density-wave interpretation that 
are often used for finding spiral arm pattern speeds. The spiral pattern is winding on a characteristic 
timcscale of ~ 500 Myrs. 

Subject headings: galaxies: fundamental parameters - galaxies: individual (NGC 1365) - galaxies: 
kinematics and dynamics - galaxies: spiral - methods: data analysis - methods: 
statistical 



1. INTRODUCTION 

The purpose of this paper is to develop and apply a 
new method of solving the Tremaine- Weinberg equations 
(Tremaine & Weinberg 1984, hereafter TW84) for a pat- 
tern speed, Q p , that is allowed to vary with radius. The 
pattern speed is an important astrophysical parameter 
for understanding the role of spiral and bar structures in 
the formation, dynamics, and evolution of a galaxy. Its 
radial behavior is not well understood. 

The method is capable of checking fundamental as- 
sumptions in the density-wave interpretation of spiral 
structure that are often used for measuring f2 p . This 
interpretation assumes that J7 p is constant (i.e., indepen- 
dent of radius ) , and that resonant radii are identifiable 
on a rotation curve of the galaxy (see Bertin & Lin 1996; 
Bertin 2000). Resonant radii are assumed to occur when 
Sip = fi, the material speed, and when Q p = SI ± n/m, 
where n is the epicycle frequency of a stellar orbit and 
m is the mode of the wave corresponding to the number 
of arms l (|B. Lindblad 1938ft . These resonances are com- 
monly referred to as corotation and Lindblad resonances 
respectively, and are identified by interpreting photomet- 
ric and kinematic features of a galaxy. The assumption 
that rip is constant in the spiral arms has yet to be con- 
firmed by direct observation. 

An alternative way of measuring Sl p is to solve the 
mo del- independent TW84 equations. These are inte- 
grated forms of the continuity equation that relate Sl p 
to the brightness of a pattern tracer and the observable 
line-of-sight velocity. When Sl p is allowed to vary with 
radius, integrating along paths that are parallel to the 
kinematic major axis of the galaxy produces an integral 



equation for Sl p , 
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where / is the specific intensity of a pattern tracer and 
v y is the line-of-sight velocity times the cosecant of the 
disk inclination. An illustration showing how a galaxy's 
disk is oriented in the coordinate system of Equation 
(1) is shown in Figure 1. In this coordinate system, y 
is measured in the plane of the disk. Equation (1) is 
similar to Equation (4) in TW84 with the exception that 
Sip is allowed to vary with radius. Although allowing 
Sl p to vary with radius may violate assumption 2. in 
TW84, it can be shown that the derivation by TW84 is 
generalizable to a pattern speed that varies with radius 
(for a proof see Engstrom 1994, hereafter E94). The 
steps involved in deriving Equation (1) for Hi are shown 
in Section 2.2. 

The new solution method uses least squares to solve 
Equation (1) for simple functional forms of Sl p . The two 
forms solved for in this paper are a polynomial in r = 
\/x 2 + y 2 , 
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and a polynomial in 1/r, 
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where on are unknown constants to be determined in 
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Fig. 1. — Illustration showing how a galaxy disk is oriented in the 
coordinate system of Equation (1). The origin of the coordinate 
system is at the galaxy's kinematic center. The kinematic major 
axis of the disk is along the y = axis. An example of an integration 
path at y = yj is shown as a solid line across the disk. 

the solution. Both equations can represent a constant 
Hp when n — 0. Although they are simple in their 
form, Equations (2) and (3) are adequate for investigat- 
ing whether fl p has a radial dependence. If evidence for 
such a dependence is found, then the pattern is shear- 
ing, and the solution obtained when n = is the mean 
pattern speed. 

To find the least squares solution, first combine Equa- 
tions (1) and (2) to form a sum on the left hand side of 
Equation (1), 



E 

i=0 



+ 00 Q 



dx 



dy 



I(x,y)v y (x,y)dx, (4) 



and similarly for Equation (3), 
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What remains inside the integrals are observable quanti- 
ties that can be calculated at different y = yj distances 
from the kinematic major axis (see Figure 1). For several 
independent calculations let, 
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for Equation (4), and similarly, 
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for Equation (5), where that matrix indexing is used for 
j rows and i columns. Also let, 



+00 



d_ 

dy 



I(x,y)v y (x,y) 



dx, 



(8) 



for the right hand sides of Equations (4) and (5). In 
Equations (6) and (7) the ith column corresponds to the 
power of r in O p . In Equations (6), (7), and (8) the 
jth row corresponds to a calculation at yj. Using these 
matrices, Equations (4) and (5) in matrix form are, 

Ga = d, (9) 

where that a = on are the unknown coefficients. Multi- 
plying both sides of Equation (9) by the transpose of G 
produces the normal equations, 



G J Ga = G J d, 



(10) 



for the least squares problem. The maximum likelihood 
solution for the best-fitting coefficients is obtained by 
multiplying both sides of Equation (10) by (G T G) _1 . 
Standard statistical tools can then be used to analyze 
the results and determine whether a constant £l p , or one 
that varies with r, provides a better solution. 

The rest of this paper is dedicated to applying the so- 
lution method to the spiral pattern of NGC 1365 and 
analyzing the results. In Section 2 is a description of 
NGC 1365 and the details of applying the method. In 
Section 3 is an explanation of the statistics used for es- 
timating the uncertainties and analyzing the results. In 
Section 4 are the results and their statistical analysis. In 
Section 5 the method is applied to different regions of 
NGC 1365 to rule out selection bias in the results. For 
all of the regions the method was applied to, the results 
are very similar to fl. In Section 6, therefore, the method 
is tested using galaxies without spiral patterns to see if 
their results are also similar to ft. In Section 7 other 
solution methods are applied to NGC 1365 and the re- 
sults are compared to those obtained using the current 
method. In Section 8 the results are compared with pre- 
vious estimates of Sl p . In Section 9 the relevance of the 
results to theories of spiral structure is discussed, and in 
Section 10 is a summary. 

2. APPLICATION TO NGC 1365 

2.1. Description of NGC 1365 

NGC 1365 is a grand-design spiral galaxy of type 
SB(s)tJ (|de Vaucouleurs et al. 1991. hereafter RC3D . An 
optical image of NGC 1365 is provided in Figure 2 to 
assist in its description. The image is reproduced from 
the Southern Digitized Sky Survey of the Space Tele- 
scope Science Institut^- The highly symmetric, coher- 
ent, two-arm pattern of NGC 1365 is ideal for the first 
application of the solution method. The disk inclina- 
tion is small enough to observe the spiral structure, yet 
large enough to accurately measure v y . The assumed 
position, orientation, and systemic velocity, V sys , that 

Original plate material is copyrighted by the Royal Observa- 
tory Edinburgh and the Anglo-Australian Observatory. The plates 
were processed into the present compressed digital form with their 
permission. The Digitized Sky Survey was produced at the Space 
Telescope Science Institute under US Government grant NAG W- 
2166. 



11.1 



The Shearing Spiral Arms of NGC 1365 



3 




TABLE 1 

Adopted Parameters of NGC 1365 



Fig. 2. — Optical image of NGC 1365 in the blue part of the 
optical spectrum. The original image is a 60 minute Illa-J exposure 
using the UK 48-inch Schmidt Telescope and a 4680 A GG 395 
filter. North is up and East is to the left. The receding half of the 
galaxy is to the bottom right. 



were used for applying the method are shown in Ta- 
ble 1. These were adopted from the analysis of the H 
I velocity field by Jorsater & van Moorsel (1995, here- 
after JvM95). NGC 1365 be longs to the Fornax clustej 
(|Burbidge fc Burbidge I9 60), but has no known close 
companions. This archetypical barred spiral galaxy is 
well studied. Fo r more informat i on ab out NGC 1365, 
see the review by IP. O. Lindbladl (fl99l . 

2.2. Hi as a Pattern Tracer 

TW84 derive pattern speed equations from a form of 
the continuity equation that assumes mass conservation. 
This assumption is difficult to satisfy for tracers of spiral 
patterns; starlight can be obscured by dust, and the ISM 
is multi-phased. When using a component of the ISM 
such as Hi, an argument is usually made for dominance 
in the abundance of the tracer within dynamical (orbital) 
timescales so that conversion among the different phases 
of the ISM is negligible. These arguments were made for 
Hi by Bureau et al. (1999) and Meidt et al. (2009, here- 
after M09), and for H 2 by Zimmer et al. (2004), Rand 
k Wallin (2004), Merrifield et al. (2006, hereafter M06), 
Meidt et al. (2008b), and Meidt et al. (2009). For the 
latter group, a known conversion factor from CO surface 
brightness to H 2 abundance is also assumed. M09 also 
combined measurements of H I with estimates of H 2 in an 
attempt to account for the total abundance of gas. The 
TW84 equations have also been applied to H-a emission 
by Hernandez et al. (2005), Emsellem et al. (2006), Fathi 
et al. (2007), Chemin & Hernandez (2009), Gabbasov et 
al. (2009), and Fathi et al (2009), the first and the last 
of which acknowledge that ionized gas neither dominates 
the ISM, nor persists within dynamical timescales, but 
claim that it may come close enough. 

The TW84 equations are derivable, however, from the 



Kinematic center R. A. (J2000) 
Kinematic center DEC. (J2000) 
Inclination angle 
Position angle 



03^ 33 m 36.4 s ± 4" 
-36° 08' 25.5" ± 4" 
40° ± 2° 
220° ± 2° 
1632 ± 3 km s" 1 



general form of the continuity equation, 

d T , . d T , . , . 
—I(x,y,t) + —I{x,y,t) v x (x,y,t) 

d 

+ dy I(yX ' V " ^ Vy ^ X ' Vl l ' = S ^' y ' ^' ^ 

that includes a source function, S, to account for the 
creation and destruction of material. This was shown by 
Westpfahl (1998, hereafter W98), who assumed that the 
intensity has a solution of the form, 

I(x,y,t)=I(x,y,t)+ f S(x,y,t')dt\ (12) 
Jo 

where the time dependance of / is due entirely to rota- 
tion. The integral of the source function is differentiable 
using Leibnitz's rule, 

J^jf S(x,y,t')dt' = S(x,y,t) + J ^S(x,y,t')dt\ 



(13) 



and from TW84 the time derivative of J is, 

^I(x, y, t) = n p (r) {v^i{x, y, t) - x^-i 



(14) 



Combine Equations (11) through (14), and the result, 

ttp(r)[y-^I(x,y,t) - x—I(x,y,t) j 
d d 

rt q 



= -J Q ls(x,y,t') dt ' 
= S(x,y,0)-S(x,y,t),(15) 



is different from that derived in TW84 by the change in 
the source function on the right-hand side. For the pur- 
pose of this paper, it is assumed that the source function 
is approximately constant within dynamical timescales. 
The result, then, is the same pattern speed equation de- 
rived by TW84, but from a form of the continuity equa- 
tion that does not assume mass conservation. With the 
change in the source function set to zero, integration in 
x removes the velocity term in that direction due to the 
boundedness of /, and produces Equation (1). 

Of the spiral pattern tracers that are available, Hi has 
several advantages: 

1.) Radio interferometers such as the Very Large Ar- 
ray (VLA) of the National Radio Astronomy Observa- 
tory (NRAO) are capable of producing high-resolution 
specific intensity and velocity maps of Hi from the 21- 
cm line. 
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2. ) To a good approximation, the 21-cm line is opti- 
cally thin, unlike, for example, the 2.6-mm line of CO. 
Assuming the 21-cm line is optically thin, it is straight- 
forward to measure the total abundance of Hi, which is 
necessary for applying the continuity equation. 

3. ) The spiral patterns of galaxies often extend much 
farther in the 21-cm line than they do in other tracers 
such as optical starlight, H-a emission, or the 2.6-mm 
line. 

4. ) Unlike optical starlight or H-a emission, the 21-cm 
line is not vulnerable to obscuration and reflection by 
dust. 

5. ) Bar patterns are often undetected in Hi maps of 
barred spiral galaxies, allowing for a separate measure- 
ment of the spiral arm pattern speed. 



2.3. Data 

Specific intensity and velocity maps made by JvM95 
were used to make maps of I and v y as described in the 
Introduction. They mapped the Hi to high resolution 
from VLA observations in the BnA, CnB, C, and DnC 
configurations, amounting to a total of 53.8 hours on 
source. The major and minor axis at FWHM of the 
synthesized beam was 11.55" x 6.32" with a position 
angle of 5.33°. The RMS noise in a line- free channel was 
0.27 mJy beam -1 . The channel width was 20.84 km s -1 . 
For more information about the observations and data 
reduction see JvM95. 

The original maps were modified using the Astro- 
nomical Image Processing System (AIPS) developed and 
maintained by the NRAO. This included using the AIPS 
task OGEOM for image translation and rotation, and 
the AIPS task MATHS for subtracting V sys from the ve- 
locity map and then multiplying it by the cosecant of the 
disk inclination. Figure 3 shows the centered and rotated 
specific intensity map overlaid with contours of constant 
v y . The pixel increment spans 2" in the x direction and 
2.61" in the y direction. Estimates of the uncertainty per 
pixel in the maps of / and v y are explained in Section 
3.1. They are 07 = 0.54 mJy km s -1 pixel -1 , and a Vy — 
12.86 km s -1 pixel -1 , respectively. 



2.4. Implementation 

The solution method was implemented within MAT- 
LAE0. Calculations of Equations (4) and (5) were re- 
stricted to the region \y\ ^ 400". Different yj were spaced 
5 pixels apart, which amounts to a separation distance 
that is 1.49" larger than the FWHM of the synthesized 
beam in that direction. This provided 62 independent 
calculations of Equations (4) and (5). To check for bias 
from this choice of spacing, larger spacings of 10 pixels 
were tried and consistent results were obtained. Differen- 
tiation was performed us ing the Savitzky-Golay method 
(|Savitzkv fc Golav 1964T ) . Integration was performed by 
summing over a row of pixels in a map of the integrand. 
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Fig. 3. — Rotated H I specific intensity map of NGC 1365 overlaid 
with contours of constant v y . The galaxy is oriented in the same 
way as shown in Figure 1. The central, most vertical contour is v y 
= km s _1 . Velocity contours are spaced 80 km s — 1 apart. The 
receding half of the galaxy is to the right. The peak flux is 250 Jy m 
s — 1 beam -1 . In the bottom left corner is a plot of the synthesized 
beam. Note that the HI does not trace the bar pattern. 



3. STATISTICS 
3.1. Uncertainties 

The analysis of the results rely on estimates of the un- 
certainties for the coefficients and the total uncertainties 
for calculations of Equations (4) and (5). The latter are 
not shown explicitly in this paper, but are part of the chi- 
squared (x 2 ; and reduced chi-squared, xt) values used in 
the analysis. How these uncertainties were estimated is 
explained below. 

The uncertainties for the coefficients are reported as 
the half-width, HW, of their 95% confidence intervals. 
The 95% confidence interval for the ith coefficient is, 



i t/ (2.5%)SE i sc cti sc £„(97.5%)SEi 



(16) 



where t v is the percentile of a Student's t - distribution 
with v degrees of freedom, and SE^ is the standard error 
for the ith coefficient (see Ramsey & Shafer 2002). The 
solutions in this paper have enough degrees of freedom 
that t„(2.5%) = £497.5%). The HW of a 95% confidence 
interval for the zth coefficient is, 

HW, = ^(97.5%)SE l . (17) 

In a similar fashion, the HW of a 95% confidence band 
for Q p when Equations (2) or (3) are summed to order n 
^ 1 is, 



HW(r) = ^(97.5%)SE(r), 



(18) 



2 The M athworks. Version 
http://www.mathworks.com 



7.10.0.499 (R2010a), 



where SE(r) is the standard error of Q p as a function of 
radius. For the solutions in this paper t„(97.5%) « 2. 

The standard errors were found using Monte Carlo 
methods. Monte Carlo methods are preferred to those 
using the covariance matrix because such methods as- 
sume <tq. i = 0, which is not true in this case. For ten- 
thousand iterations, normally distributed random noise 
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with standard deviations of gq . 4 and were added 
to G and d respectively, and solutions were found for 
the coefficients. The standard deviations in the distri- 
butions of the results were adopted as estimates of SEj. 
Likewise, the standard deviations for distributions of fi p 
over a range of r were adopted as estimates of SE(r). 
Monte Carlo methods produced confidence intervals and 
confidence bands ~ 1 - 10% greater than methods using 
the covariance matrix. The difference was greatest for 
forms of fl p having the most terms. 

Estimates of aQ j i and a^. were found by propagating 
estimates of the uncertainty per pixel, 07 and a Vy given 
in Section 2.3, as well as the uncertainties in x and y, 
through the calculations of G and d. The RMS noise in a 
line- free channel is a lower limit to 07 in a map of / where 
emission is present, so doubling the RMS noise should 
adequately account for tr/. To estimate a Vy , the uncer- 
tainty per pixel in velocity resolution, the uncertainty in 
V sys , and the uncertainty in inclination were propagated 
throu g h the calculation of v y . Following the reasoning of 
iRots I (|1974f ) , the uncertainty per pixel that is due to the 
velocity resolution was estimated as the channel width 
divided by the number of channels in which a detection 
occurs. JvM95 show that only the brightest Hi appears 
in more than two channels of the data cube, the weakest 
detections appear in one, and the majority appears in 
about two. The uncertainty due to the velocity resolu- 
tion was therefore set to one-half the channel width. The 
uncertainties in x and y were calculated from the uncer- 
tainties in the kinematic center and the position angle 
(PA) . The uncertainty in y also included the uncertainty 
in the inclination angle. 

To avoid the need for propagating the uncertainties in 
x and y that are due to uncertainties in the distance to 
NGC 1365, the instrumental units of the VLA, [ftp] = km 
s _1 arcsec -1 , are used for calculating Equations (4) and 
(5) and reporting the results. Furthermore, using units 
common for galactic disks, [Q p ] = km s _1 kpc -1 , would 
introduce an unnecessary amount of uncertainty when 
determining the functional form of fl p that provides the 
best solution, or when comparing the results with those 
obtained using other solution methods and from previ- 
ous estimates. After the analysis, the most important 
results are reported in units commonly used for galac- 
tic disks. The uncertainties for the converted results are 
reported as la uncertainties, which were estimated by 
propagating the standard errors for the coefficients and 
the uncertainty for the distance through the conversion. 

Errors in the PA can contribute to errors in the coeffi- 
cients that are not represented in the propagated uncer- 
tainties. Monte-Carlo error estimations were used to ex- 
amine this effect. For ten-thousand iterations, a PA was 
randomly chosen from a normal distribution of angles 
centered on the assumed PA whose standard deviation 
was equal to the uncertainty in PA, apA- The result- 
ing distributions in the coefficients were typically non- 
gaussian, and produced pattern speeds that were most 
often smaller than those derived using the assumed PA. 
It would therefore be inappropriate to add the standard 
deviations in the results for the Monte-Carlo estimates 
in quadrature to the uncertainties in the coefficients in 
order to account for a pa- Instead, the effect of an in- 
correct PA is discussed for each galaxy on a case by case 



basis. The percent error due to an incorrect PA is much 
smaller than those that are typically reported for the 
original solution method of TW84. The percent error for 
the original solution method can be as large as 25% for a 
2° error in PA, and 100% for a 5° error in PA (Debattista 
2003; Debattista & Williams 2003). The precent errors 
for the current method are about 1/4 as large. 

Finally, the total uncertainty for each jth calculation 
of Equation (4) or (5), aj, was used for finding y 2 ( an d 
xt)- These were estimated as, 

n 

by following the reasoning of Bevington & Robinson 
(2002) to take into account both and OQ j A . 

3.2. Analysis 

The primary goal of the analysis is to determine the 
functional form of fl p that provides the best solution, 
and thus the best estimate of Q p . Starting at order n 
= 0, the value of xt was compared with the value ob- 
tained after the addition of an extra term. If there was 
no improvement in xti then the additional term did not 
provide an improvement to the solution. If there was 
an improvement in xt, then the F test of an additional 
term was used to find the statistical significance of the 
improvement. It consists of calculating 

Ay 2 

F = ^-, (20) 

xi 

and finding the probability, P, of F being as large with 
the addition of an extra term, should the coefficient of 
the extra term be zero (Bevington & Robinson 2002). A 
value of P < 1% is convincing evidence that the extra 
term is non-zero, a value of P ~ 5% is suggestive, but 
inconclusive, and larger values are an indication that the 
improvement is statistically insignificant. A lower limit 
of 0.01% was chosen as a cutoff for reporting P. Such 
a small percentage means that the improvement in xt 
is extremely statistically significant. More information 
about h ow to interpret values of P in the F test can be 
found in lRamsev fc Schafer I (|2002[ ). For the n = 1 term, 
an F test providing convincing evidence that it is non- 
zero is equivalent to relaxing the assumption that Q p is 
constant. 

Terms were added based on the results of the F test 
until there was no longer an improvement in \t > the im- 
provement was no longer statistically significant, or G T G 
became ill conditioned for numerical inversion. The ma- 
trix G T G will eventually become ill conditioned for nu- 
merical inversion in the same manner that is expected 
for least squares solutions using polynomials with suc- 
cessively higher order terms. When this occurs the solu- 
tion is unstable and t he overall statistics are unreliable 
(|see Aster et al. 2005T) . 

In addition to the above statistics, R 2 was found for a 
solution when Equation (2) or (3) was summed to order 
n ^ 1. In this paper R is the percent of the variance for 
the solution when n = that is explained by the solution 
using more than 1 term. The value of R 2 in combination 
with the results from the F test were used for comparing 
the goodness of fit for the two different forms of il p . 
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TABLE 2 

Best-Fit Coefficients for NGC 1365 







a ± HW 


qi ± HWi 


a 2 ± HW 2 


a 3 ± HW 3 


Equation 


n 


(km s _1 arcscc -1 ) 


10~ 2 (km s _1 arcscc~ 2 ) 


10~ 4 (km s _1 arcscc~ 3 ) 


10~ 7 (km s _1 arcsec -4 ) 


(2), (3) 





1.27 ± 0.10 








(2) 


1 


2.68 ± 0.38 


-0.55 ± 0.13 






2 


4.05 ± 0.83 


-1.65 ± 0.57 


0.20 ± 0.09 






3 


5.93 ± 1.43 


-4.06 ± 1.64 


1.11 ± 0.59 


-1.05 ± 0.64 






10~ 1 (km s _1 arcscc -1 ) 


10 2 (km s- 1 ) 


10 3 (km s — 1 arcsec) 




(3) 


1 


-0.59 ± 2.88 


2.87 ± 0.61 






2 


-1.04 ± 4.28 


3.03 ± 1.29 


0.95 ± 6.65 





4. RESULTS AND ANALYSIS 

In this section the results and analysis are presented 
by the equation used for O p , and then using plots. At 
the end of this section conclusions are drawn from the 
analysis. The results for the best-fit coefficients, ai, and 
the half- widths of their 95% confidence intervals, HWi, 
are shown in Table 2. The statistics for analyzing the 
results are shown in Table 3. 

4.1. Results From Equation (2) 

Solutions using Equation (2) show improvements in xl 
that are statistically significant to order n — 3 accord- 
ing to the F test (P < 0.01%). The overall statistics, 
however, are only reliable to order n = 2, after which 
point G T G became ill conditioned for numerical inver- 
sion. The improvement in xl for order n — 1 is extremely 
statistically significant (P < 0.01%), and explains more 
than half of the variance (R 2 — 59.5%) for the solution 
when n = 0. The improvement in xl f° r order n = 2 
is as statistically significant (P < 0.01%), and explains 
more of the variance (i? 2 = 70.1%) for the solution when 
n = 0. Since there is convincing evidence for a radial 
dependence in f2 p , the solution when n = is the mean 
pattern speed. It follows that the i — 1 coefficient for 
the n = 1 solution is the mean shear rate, and the i = 2 
coefficient for the n = 2 solution is a correction for the 
radial variance in the shear rate. 

4.2. Results From Equation (3) 

The solution using Equation (3) to order n = 1 shows 
an improvement in xl that is extremely statistically sig- 
nificant according to the F test (P < 0.01%). Its value of 
F = 525 is quite large, and is more than twice the value 
of F — 241 for the n = 1 solution using Equation (2). 
This solution also explains more of the variance (R 2 = 
79%) for the solution when n — than do the solutions 
using Equation (2). The i = term for this solution is 
indistinguishable from zero given the size of its uncer- 
tainty. To a very good approximation, the i — term 
is negligible, and the form of this solution is 1/r. The 
n = 2 solution did not provide a statistically significant 
improvement in xl (P = 7.57%). 

4.3. Plots of the Results 

Plots of il p from the results in Table 2 are more reveal- 
ing. Figure 4 shows plots of adopted from the rotation 
curve of JvM95, possible locations for 2 and 4 arm Lind- 
blad resonances, and fi p from the results for the solution 
using Equation (2) to order n = 2. The purpose of this 



TABLE 3 
Statistics for NGC 1365 











P 


R 2 


Equation 


n 


xl 


F 


(%) 


(%) 


(2), (3) 





21.9 








(2) 


1 


4.43 


241 


< 0.01 


59.5 


2 


2.30 


56.7 


< 0.01 


70.1 




3 


0.72 


129 


< 0.01 


76.3 


(3) 


1 


2.28 


525 


< 0.01 


79.0 


2 


2.19 


3.37 


7.14 


79.0 



figure is to build a complex plot one layer at a time, the 
format of which is used in subsequent plots of the results. 
The uncertainties for are omitted because they are as 
small or smaller than the circles representing it. Results 
are not plotted for r < 50", wherein there is not an Hi 
pattern. 

Figure 5 shows plots for all of the results in Table 2. 
Panel (a) shows the result when n = 0. Panels (b), (c), 
and (d) show pattern speeds having the form of Equation 
(2) to order n = 1, 2, and 3 respectively. Panels (e) and 
(f) show pattern speeds having the form of Equation (3) 
to order n = 1 and 2 respectively. Upon comparison with 
the other solutions, the solution in panel (a) is clearly 
the mean pattern speed. The larger confidence bands 
in panel (d) for the n = 3 solution using Equation (2) 
are a result of the instability of that solution. There 
is hardly any noticeable difference between panels (e) 
and (f). Note the progression in panels (a) thorough (d) 
where the form of Q p converges to a form resembling that 
in panels (e) and (f). As terms are added to Equation 

(2) , the radial behavior of f2 p becomes increasingly more 
concave up. The same behavior is achieved by Equation 

(3) with fewer terms. 

Plots of the results do not show the pattern speed of the 
bar. There are enough terms in the forms of il p in panels 
(d) and (f) to produce a turnover that is concave down 
around 100"-150", inside of which would be a constant 
flp for the bar. Although NGC 1365 is often regarded as 
a typical case for a continuous bar and spiral pattern hav- 
ing the same pattern speed, the bar pattern speed may be 
decoupled from the spiral arm pattern speed (Sellwood 
& Sparke 1988, Sellwood 1993 and references therein). 

To explore the effect of an incorrect PA, solutions were 
found for a range of PA within ± 2a pa of the assumed 
PA. The results for the mean pattern speed and the val- 
ues of F for the n = 1 solutions are shown in Figure 6. 
The changes in the mean pattern speed are indicative of 
the changes for all forms of £l p in Figure 5. When the 
mean pattern speed decreases or increases by a certain 
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Fig. 4. — Plots for explaining how subsequent plots of the results are formatted. Panel (a) shows the material speed, f2, represented by 
circles. Panel (b) shows S7 with possible locations for Lindblad resonance. Dotted and dash-dot lines represent possible locations for 2 and 
4 arm Lindblad resonances respectively. Panel (c) shows the pattern speed, Q p , using the results for Equation (2) to order n = 2 and its 
95% confidence bands. These are represented by solid and dashed lines respectively. Panel (d) is a combination of panels (b) and (c). See 
text for more details. 



amount, the other forms generally do so as well. If the 
assumed PA is incorrect, the results are most likely over 
estimating il p . The difference could be as large as ~ 0.2 
km s _1 arcsec -1 (or ~ 2.3 km s -1 kpc -1 using the as- 
sumed distance stated in Section 4.4). The mean percent 
difference within ± lap a is 7%. The largest percent dif- 
ference within ± Iupa is ~ 18%. For all of the changes 
in PA shown, there is convincing evidence for shear in 
the pattern, and the values of F are the largest for O p 
having form of Equation (3). 

4.4. Conclusions 

The best estimate of VL P is adopted from the results for 
Equation (3) to order n = 1. This conclusion is based 
on three features of the results: 1) this solution has the 
greatest values of F and R 2 ; 2) the results for Equation 
(2) converge towards a form resembling this solution as 
successive terms are added; and 3) this solution has a 
simpler form than the n — 2 solution using Equation 
(2). Although inferring the exact functional form of O p 
is beyond the scope of the results, its radial behavior is 
approximately 1/r. 

There are no clear indications of resonances associated 
with the spiral arms. Panel (e) in Figure 5 does not 
show a single point of corotation resonance. In the same 
panel, Lindblad resonances can neither be confirmed nor 
ruled out beyond 350". The intersection of f2 p and the 
possible locations for Lindblad resonances beyond 350" 



TABLE 4 

Results for NGC 1365 in Common Units 



Mean pattern speed 
Mean shear rate 
Best estimate of Q v 



14.5 (± 2.16) km s~ 
- 0.71 (± 0.24) km s" 
2.87 (± 0.61) x 10 2 km s" 1 jr 



kpc 
1 kpc" 



may be a coincidence due to the precision of the results. 
The results exclude nonlinear mode coupling of the bar 
and spiral arms, which require inner Lindblad resonances 
where the spiral arms begin (Tagger at al. 1997; Masset 
& Tagger 1997). The effect of an incorrect PA would 
not increase fi p enough to account for the lack of inner 
Lindblad resonances at the beginnings of the spiral arms. 

The most important results are summarized in Ta- 
ble 4 using common units. The units were converted 
using an assumed distance of 18.1 (± 2.6) Mpc to 
NGC 1365, which was estimated from the mean of 37 
redshift-independent measurements found in the Nasa 
Extragalactic Database (NED). The la uncertainties are 
shown in parenthesis. About half of the la uncertainties 
are due to the uncertainty in the distance. 

5. A TEST FOR SELECTION BIAS 

There is a possibility that the results for NGC 1365 
are biased by selecting the region \y\ < 400" (hereafter 
region A) for applying the method. To test for whether 
this occurred, additional solutions were found for the re- 
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Fig. 5. — Plots of pattern speeds from the results in Table 2. They are formatted in the same way as panel (d) in Figure 4. See text for 
details. Note the similarity of panels (d), (e), and (f), with larger confidence bands in panel (d). The material speed fits within the 95% 
confidence bands in panels (e) and (f). 



gions \y\ < 200", -400" sC y ^ 0", and 0" < y < 400" 
(hereafter regions B, C, and V respectively). Solutions 
for region B use less of the Hi in the outer part of the 
disk than solutions for region A. Comparing the results 
for these two regions will show whether Q p behaves much 
differently in the inner and outer parts of the disk. Com- 
paring the results for regions C and V will show whether 
il p behaves much differently for either side of the galaxy. 
The results for regions B, C, and V are shown in Table 5 
up to the same order as the results for region A in Table 
2. 

The coefficients for regions £>, C, and T> are consistent 
with those for region A. This is easily seen by comparing 
plots of f2 p , which are shown in Figures 7, 8, and 9 for 
regions £>, C, and V respectively. They show the same 
progression as was observed in Figure 5 where that the 
form of ftp in panels (a) through (d) converges to a form 



resembling that in panels (e) and (f). Any differences 
among Figures 5, 7, 8, and 9 are very small. Compared 
with Figure 5, the plotted pattern speeds in Figures 7, 8, 
and 9 have slightly larger confidence bands. This is due 
to having fewer calculations of Equations (4) and (5). 

To see whether the small differences in the plots are 
significant, the pattern speeds in panel (e) of the figures 
were subtracted from each other and plotted. Figure 10 
shows four cases to consider: in panel (a) is the result 
for region B subtracted from the result for region A; in 
panel (b) is the result for region C subtracted from the 
result for region A; in panel (c) is the result for region T> 
subtracted from the result for region A; and in panel (d) 
is the result for region V subtracted from the result for 
region C. The half- widths of the confidence bands shown 
in the plots are la uncertainties that were estimated by 
adding the standard errors in quadrature. The largest 
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TABLE 5 

Best-Fit Coefficients for Different Regions of NGC 1365 



Equation 


n 


Region 


a ± HW 

(km s -1 arcscc -1 ) 


ai ± HWi 
10 -2 (km s -1 arcscc -2 ) 


a 2 ± HW 2 
10~ 4 (km s -1 arcscc -3 ) 


a 3 ± HW 3 
10 -7 (km s -1 arcsec -4 ) 


(2), (3) 





B 
C 
V 


1.37 ± 0.12 
1.20 ± 0.18 
1.31 ± 0.13 








(2) 


1 


B 
C 
V 


2.75 ± 0.50 
3.23 ± 0.66 
2.39 ± 0.45 


-0.57 ± 0.20 
-0.74 ± 0.23 
-0.44 ± 0.16 








2 


B 
C 
V 


4.65 ± 1.20 
4.42 ± 1.38 
3.81 ± 1.44 


-2.27 ± 1.02 
-1.86 ± 1.17 
-1.50 ± 0.97 


0.34 ± 0.20 
0.23 ± 0.24 
0.17 ± 0.14 






3 


B 
C 
V 


7.12 ± 2.10 
5.37 ± 2.18 
6.35 ± 2.28 


-5.95 ± 2.86 
-3.25 ± 2.93 
-4.62 ± 2.39 


1.99 ± 1.24 
0.83 ± 1.25 
1.31 ± 0.81 


-2.28 ± 1.74 
-0.81 ± 1.70 
-1.27 ± 0.87 








10 — 1 (km s -1 arcscc -1 ) 


10 2 (km s- 1 ) 


10 3 (km s _1 arcsec) 




(3) 


1 


B 
C 
V 


0.40 ± 3.66 
-0.06 ± 3.83 
-1.39 ± 4.53 


2.70 ± 0.71 
2.78 ± 0.78 
3.03 ± 0.97 








2 


B 
C 
V 


0.55 ± 5.34 
-1.25 ± 6.17 
-0.74 ± 6.41 


2.67 ± 1.50 
3.23 ± 1.98 
2.79 ± 2.03 


0.17 ± 7.18 
2.35 ± 9.28 
1.75 ± 13.4 
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Fig. 6.— Effect of an incorrect PA for NGC 1365. Panel (a) 
shows how the mean pattern speed, represented by points, varies 
with PA. For reference, the solid and dashed lines in panel (a) 
represent the mean pattern speed and its 95% confidence interval 
obtained using the assumed PA. Panel (b) shows the values of F 
for the n = 1 solution as open squares for Equation (2) and solid 
diamonds for Equation (3). The dashed line in panel (b) is the 
value of F at which P = 1%. 



differences, as well as the largest uncertainties, occur for 
smaller radii. The plots show that the differences are 
less than 10%, but given that zero falls well within the 
confidence bands, they are indistinguishable from zero, 



TABLE 6 

Statistics for Different Regions of NGC 1365 













P 


n 2 


Equation 


n 


Region 


xl 


F 


(%) 


(%) 


(2), (3) 





B 


5.58 












C 


7.59 












V 


38.6 








(2) 


1 


B 


2.39 


39.7 


< 0.01 


52.6 




C 


1.86 


93.6 


< 0.01 


79.2 






V 


3.34 


317 


< 0.01 


46.6 




2 


B 


1.38 


21.4 


< 0.01 


68.4 






C 


1.29 


13.9 


0.09 


85.4 






V 


3.21 


2.23 


14.7 


55.7 




3 


B 


0.62 


34.0 


< 0.01 


75.0 






C 


0.80 


18.0 


0.02 


87.0 






V 


0.74 


94.1 


< 0.01 


68.7 


(3) 


1 


B 


1.14 


114 


< 0.01 


76.5 




C 


1.97 


86.5 


< 0.01 


86.0 






V 


2.22 


493 


< 0.01 


71.6 




2 


B 


1.18 






76.5 






C 


1.98 






86.5 






V 


2.55 






71.7 



and therefore insignificant. 

The statistics for the solutions using regions £>, C, and 
V are shown in Table 6 for comparison with those for 
region A shown in Table 3. Like those in Table 3, the 
statistics in Table 6 for the solutions using Equation (2) 
to order n = 3 are unreliable because at this point G T G 
became ill-conditioned for numerical inversion. All of the 
F test results for region B are consistent with those for 
region A. The outer disk is not biasing the analysis of the 
results. The F test results for regions A, C, and T> are 
consistent with each other up to order n = 1 for Equation 
(2), and up to order n = 2 for Equation (3). For region 
C, the n = 2 solution for Equation (2) did not provide a 
statistically significant improvement in \1 (P — 6.78%). 
For region T> the n — 2 solution for Equation (2) did 
not provide an improvement in xl ■ Considering that the 
statistics for region B are consistent with those for region 
A, these differences show that the inner disk contributes 
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Fig. 7. — Plots of pattern speeds for region B. The figure is formatted the same way as Figure 5. In panel (d) the arrow at the end of 
the upper confidence band points in the direction it continues. 



the most to the concave up behavior of f2 p observed in 
the plotted pattern speeds. Despite the noted differences, 
the radial behavior of fl p is approximately the same for 
both sides of the galaxy. 

There are two consistencies among all four regions that 
are worth emphasizing. The improvement in \t f° r so ~ 
lutions using Equations (2) and (3) to order n = 1 is 
extremely statistically significant for all four regions (P 
< 0.01%). Furthermore, the values of F and R 2 for so- 
lutions using Equation (3) to order n = 1 are the largest 
among the results that provided a statistically signifi- 
cant improvement in xt- The existence of shear in the 
pattern, or the approximate 1/r behavior of fi p , is not 
biased by selecting region A for applying the method. 



6. A TEST TO VERIFY THAT THE METHOD IS 
NOT MEASURING Q 

The best estimate of fl p is very similar to ft. This raises 
the question of whether the method is actually measur- 
ing Q. To explore this possibility, a test was performed 
by applying the solution method to Holmbcrg II and IC 
2574. Both are dwarf irregular galaxies (RC3) that have 
no coherent or well organized patterns, and show mate- 
rial clumped around expanding H I bubbles (Puche et al. 
1992; Walter & Brinks 1999). Results for their pattern 
speeds would represent the averaged angular rotation of 
the individual clumps of material, and are therefore ex- 
pected to differ somewhat from Q. 

Data maps for both galaxies were obtained from the 
THINGS H I galaxy survey. Information about the obser- 
vations, dat a reduction, and para meters of the maps can 
be found in iWalter et al. I ([20081 ) . The rotation curves, 
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Fig. 8. — Plots of pattern speeds for region C. The figure is formatted the same way as Figure 5. In panel (d) the arrow at the end of 
the upper confidence band points in the direction it continues. 



as well as the position and orienta tion parameters, 
were adopted from those determined bv lS. -H. Oh et al71 
(|2010. submitted! ) for Holmberg II and by S. -H. Oh et 
al. (2008, 2010, submitted) for IC 2574. The position 
and orientation parameters are shown in Table 7. The 
rotation curves, as well as the uncertainties for the pa- 
rameters shown in Table 7, were provided by S. -H. Oh 
(2010, private communication). In the rest of this section 
the descriptions of the data, the results, and the analysis 
are presented by galaxy. These subsections are followed 
by conclusions of the test. 

6.1. Holmberg II 

Figure 11 shows the centered and rotated specific in- 
tensity map of Holmberg II overlaid with contours of 
constant v y . In the figure one can see the clumpy dis- 
tribution of the Hi, and the lack of a coherent or well 
organized pattern. Calculations of Equations (4) and (5) 



were restricted to the region \y\ ^ 400". Different yj 
were spaced 9 pixels apart, which amounts to a sepa- 
ration distance that is 0.37" larger than the FWHM of 
the synthesized beam in that direction. This provided 
36 independent calculations of Equations (4) and (5). 
Consistent results were obtained using different regions 
of implementation and larger spacings between different 

Vr 

The statistics for the results are shown in Table 8. For 
both forms of f2 p , the improvement in xt is extremely 
statistically significant to order n = 1 (P < 0.01%). Ad- 
ditional terms did not provide an improvement in xt- 
The n = 1 solution using Equation (2) has the largest 
value of F and R 2 , and is therefore adopted as the best 
estimate of H p . For this solution ao = 0.29 (± 0.12) km 
s -1 arcsec -1 and a-y = — 3.16 (± 3.76) x 10 -4 km s -1 
arcsec -2 . 
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Fig. 9. — Plots of pattern speeds for region T>. The figure is formatted the same way as Figure 5. 



Figure 12 shows plots of fl p for the solutions whose 
statistics are shown in Table 8. The best estimate of O p 
is shown in panel (b). Compared with the sloping line 
in panel (b), the result in panel (a) is clearly the mean 
pattern speed. There are more negative residuals in all 
of the panels. All of the panels also show the pattern 
corotating with the material in the inner part of the disk, 
and rotating faster than the material in the outer disk. 
There may be a faster rotating clump of material in the 
outer disk that is leveraging the solution. If the method 
were measuring f2, then the result in panel (a) would be 
smaller, and the slope of the line in panel (b) would be 
larger. 

The effect of an incorrect PA was explored in the same 
way as for NGC 1365. The results are shown in Figure 
13. Most of the results for the mean pattern speed are 
smaller than the value obtained using the assumed PA. 
The mean percent difference within ± \gpa is 13%. The 



largest percent difference within ± lapA is ~ 48%. The 
pattern speed for the solutions to order n = 1 would be 
coincident with for a percent difference ~ 30%, but in 
panel (b) one can see that these solutions generally do 
not provide a statistically significant improvement to \t ■ 

6.2. IC 2574 

Figure 14 shows the centered and rotated specific in- 
tensity map of IC 2574 overlaid with contours of constant 
v y . Like Figure 11, one can see the clumpy distribution 
of the Hi and the lack of a coherent or well organized 
pattern. Calculations of Equations (4) and (5) were re- 
stricted to the region \y\ < 600". Different yj were spaced 
9 pixels apart, which amounts to a separation distance 
that is 1.11" larger than the FWHM of the synthesized 
beam in that direction. This provided 52 independent 
calculations of Equations (4) and (5). Consistent results 
were obtained using different regions of implementation 
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FlG. 10. — Plots for comparing f! p among regions A., B, C, and T>, where that Q p has the form of Equation (3) to order n = 1. The solid 
line is the difference between two regions. The dashed line is the lcr uncertainty for the difference. The dotted line at zero is provided for 
reference. See text for details. 



TABLE 7 

Adopted Parameters of Holmberg II and IC 2574 



Parameter Holmberg II IC 2574 

Kinematic center R. A. (J2000) 08 h 34 m 06.5 s ± 8.25" 10 h 28 m 27.7 s ± 7.5" 

Kinematic center DEC. (J2000) +70° 43' 24" ± 12.45" +68° 24' 59" ± 15" 

Inclination angle 49° ± 6° 55° ± 6° 

Position angle 275° ± 10° 53° ± 2° 

V SV3 156.5 ± 1 km s" 1 53 ± 1 km s" 1 



TABLE 8 
Statistics for Holmberg II 











P 


R 2 


Equation 


n 


xl 


F 


(%) 


(%) 


(2), (3) 





2.85 








(2) 


1 


1.84 


20.1 


< 0.01 


12.7 


2 


1.87 






12.7 


(3) 


1 


2.03 


15.2 


< 0.01 


9.8 




2 


2.09 






10.0 






TABLE 9 






Statistics for IC 2574 












p 


R 2 


Equation 


n 




F 


(%) 


(%) 


(2), (3) 





1.83 








(2) 


1 


2.82 






7.0 


(3) 


1 


1.92 






1.2 



and larger spacings between different yj. 

The statistics for the results are shown in Table 9. For 
both forms of il p , the n = 1 term did not provide an im- 
provement in xt- The solution when n = is therefore 
adopted as the best estimate of Q p . For this solution a 



= 0.15 (± 0.05) km s -1 arcsec -1 . The lack of convincing 
evidence for a radial dependence in tt p does not require 
that the pattern is rigidly rotating. The solution method, 
or the simple functional forms used for Cl p , may be in- 
adequate for detecting variations in the pattern speed of 
IC 2574. 

Figure 15 shows plots of tt p for the solutions whose 
statistics are shown in Table 9. Like the results for Holm- 
berg II, none of the pattern speeds shown resembles 0. 
There is a clear downward trend in fi that is not seen in 
any of the plots of Q p . Panel (a) shows the best estimate 
of £l p , which is less than ft for most of the galaxy. If 
the method is measuring f2, then the result in panel (a) 
would be larger. All three panels show the pattern coro- 
tating with the material in the outer part of the disk, 
specifically around 500". If the method is measuring ft, 
then the result in panel (a) would be larger. 

The effect of an incorrect PA was explored in the same 
way as for NGC 1365. The results are shown in Figure 16. 
All of the results for the mean pattern speed fit within the 
95% confidence interval for the value obtained using the 
assumed PA. The mean percent difference within ± 1<jpa 
is 2%. The largest percent difference within ± 1<jpa is 
~ 22%. For all of ± lcr pa, and most of ± 2a pa, the n 
= 1 solutions did not provide a statistically significant 
improvement in xt- The largest pattern speed shown in 
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y, removes the spatial derivatives. Then, Equation (1) 
reduces to, 



x (arcsec) 

Fig. 11. — Rotated HI specific intensity map of Holmberg II 
overlaid with contours of constant v y . The galaxy is oriented in 
the same way as shown in Figure 1. The central, most vertical 
contour is v y = km s . Velocity contours are spaced 26 km s — 1 
apart. The receding half of the galaxy is to the right. The peak 
flux is 599 Jy m s beam -1 . In the bottom left corner is a plot 
of the synthesized beam. 



panel (a), ~ 0.17 km s -1 arcsec -1 , would not provide 
much of a better fit to the material speed than the value 
obtained using the assumed PA. 

6.3. Conclusions 

The results for Holmberg II and IC 2574 are two exam- 
ples of when J7 p is not very similar to il. Unlike plots of 
the results for NGC 1365, plots of the results for Holm- 
berg II and IC 2574 suggest specific regions of corotation. 
Their results are consistent with patterns that are decou- 
pled from the material. If the spiral pattern of NGC 1365 
is decoupled from the material, the difference may be too 
small to distinguish with the precision of the results. It 
is also possible that the assumed PA of NGC 1365 is in- 
correct, in which case the pattern is most likely rotating 
more slowly than the material. 

The most important results for Holmberg II and IC 
2574 are summarized in Table 10 using common units. 
The units were converted using an assumed distance of 
3.33 (± 0.99) Mpc to Holmberg II and 3.85 (± 0.68) Mpc 
to IC 2574. The distances and their uncertainties were 
estimated from the mean of 6 redshift-independent mea- 
surements for Holmberg II, and 11 such measurements 
for IC 2574, all found in the NED. 

7. COMPARISON WITH THE RESULTS FROM 
OTHER SOLUTION METHODS 

7.1. The Original Solution Method of TW84 

TW84 originally sought a solution for bar pattern 
speeds, allowing them to pull Q p outside the integral in 
Equation (1). Carrying through the integration in x for 
the second term on the left-hand side of Equation (1) re- 
moves that term because of the boundedness of /. Subse- 
quent integration from y to +oo, or similarly for negative 



51 p / I(x, y) xdx = 



I(x,y)v y (x,y)dx. (21) 



Several independent calculations at different distances 
from the kinematic major axis can be used to find 
an average of the results and measure the uncertainty. 
It is common practice to normalize Equation (21) by 

f*™I(x,y) dxl ^e.s. Merrifield fc Kuilken 1995ft . 

If the spiral arms rotate with a constant f2 p , then 
the original solution method is applicable to those pat- 
terns as well. However, there are many instances where 
that the results sugges t otherwise. Resul ts for NGC 
3031 (W98), NGC 29 15| (IBureau at al. 1999ft. NGC 519 4 
(IZimmer et al. 2004ft. NGC 10681 (IRand fc Wallin 20041) 
NGC 4321| (IHernandez et al. f 2005ft . and NGC 6943 
(|Fathi et al. 2007ft show a decreasing trend for increasing 

\y\- 

Such trends are consistent with a shearing pattern. 
Consider that according to the mean value theorem for 
integration, 

f+°° f +oc ( d d \ 

J J ttp(r){y—I(x,y)-x—I(x,y)jdxdy 



Q p (f 



I(x,yj)xdx, (22) 



where f2 p (r) is the mean pattern speed for the region 
of integration, and f is locat ed somewhere in the region 
(jGradshtevn fc Rvzhik 2007ft . The lower limit for val- 
ues of f increases with increasing \y\. Equation (22) 
shows that if Q p is decreasing with radius, then the re- 
sult from calculating Equation (21) will be £l p (f), whose 
value should show a decreasing trend with increasing \y\. 

The original solution method of TW84 was applied to 
the Hi spiral pattern of NGC 1365 to see if the same 
statistics explained in Section 3 would show evidence for 
a decreasing trend in the results plotted versus |y|. Indi- 
vidual calculations of Equation (21) were spaced 5 pixels 
apart, and covered the same region as was used for the 
current method (\y\ < 400"). The half-widths of the 95% 
confidence intervals for the coefficients, and 95% confi- 
dence bands for the trend lines fit to the results, were 
calculated in a similar way as Equations (17) and (18), 
with the exception that the standard errors were esti- 
mated from the covariance matrix. 

Figure 17 shows constant and sloping trend lines fit to 
the results. The constant trend line at 0.80 (± 0.22) km 
s -1 arcsec -1 in panel (a) is less than the mean pattern 
speed of 1.27 (± 0.10) km s -1 arcsec -1 obtained using 
the current method. The difference is due to the lower 
value of Q p in the outer disk contributing to O p (f) for all 
calculations. Note that in panel (a) there are more posi- 
tive residuals for smaller \y\ and more negative residuals 
for larger \y\, which occurs separately for both positive 
and negative y. The sloping trend line in panel (b) has 
an intercept on the vertical axis at 1.86 (± 0.30) km s -1 
arcsec -1 and a slope of — 0.53 (± 0.13) x 10 -2 km s -1 
arcsec -2 . An F test showed that the sloping trend line 
provides a statistically significant improvement to xt (P 
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Fig. 12. — Plots of pattern speeds from the solutions for Holmberg II. The figure is formatted the same way as Figure 5. Panel (a) shows 
the mean pattern speed. Panels (b), and (c) show pattern speeds having the form of Equation (2) to order n = 1, and 2 respectively. 
Panels (d) and (e) show pattern speeds having the form of Equation (3) to order n = 1 and 2 respectively. 



= 1.32%). 

7.2. Allowing for Q p to Vary With Radius 
7.2.1. The Solution Method of E94 

When £l p is allowed to vary with radius, solutions of 
the continuity equation integrated over x and y are non- 
trivial. Such solutions were first explored by E94, who 
showed that a numerical solution is possible by making 
a switch to polar coordinates. The result is a Volterra 
equation of the first kind for Q p , 

r+oo r + oo 

/ n p (r)K jt i(ri,yj)dr = I(x,y 3 ) v y (x,yj) dx, 

J r—yj J — oo 

(23) 

where, 

Ki,i = ri{l(yfi yl Vj ) H-yfi yj,yj)}, (24) 



or similarly for negative y. For an illustration showing 
the radial binning in Equation (18), see Figure 1 of Mcidt 
et al. 2008a. Equation (23) is solvable using back substi- 
tution methods, but as E94 points out, such methods are 
ill-advised because of the triangular shape of K and its 
sensitivity to symmetries in the data. This was indeed 
found to be the case by M06, whose solutions for the CO 
spiral pattern of NGC 1068 became increasingly unsta- 
ble for smaller radii. They showed examples of some of 
the results for large radial binnings in the discretization 
of Equation (23) that produced a declining trend in Sl p 
with radius. 

The method of E94 was applied to the H I spiral pat- 
tern of NGC 1365 to see how the results compared with 
those obtained using the current method, but the solu- 
tions were very unstable. Radial binnings as small as 
the synthesized beam produced results that oscillated by 
up to 4 orders of magnitude. The results for larger bin- 
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TABLE 10 

Results for Holmberg II and IC 2574 in Common Units 



Galaxy 




Result 


Holmberg II 


Mean pattern speed 
Mean shear rate 
Best estimate of f2 p 


12.64 (± 4.02) km s" 1 kpc" 1 
-1.17 (± 0.80) km s" 1 kpc~ 2 
17.7 (± 6.47) km s" 1 kpc" 1 

-1.17 (± 0.80) km s" 1 kpc" 2 r 



IC 2574 Best estimate of U p 8.10 (± 1.92) km s" 1 kpc" 1 
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Fig. 13. — Effect of an incorrect PA for Holmberg II. The figure 
is formatted the same way as Figure 6. 



nings ranging from 100" - 200" oscillated by an order of 
magnitude. The wildly oscillating behavior and the lack 
of a consistent trend make these results unreliable, and 
are therefore not reported. A single radial binning from 
50" to 400", however, produced a mean pattern speed of 
1.14 (± 0.06) km s -1 arcsec" 1 for negative y, and 1.36 
(± 0.07) km s _1 arcsec -1 for positive y, in agreement 
with the results for regions C and T> respectively in Ta- 
ble 5. The uncertainties reported for these solutions may 
be larger because they do not include an estimate of the 
error introduced by rebinning the data into cylindrical 
coordinates. The average of the two results is 1.25 (± 
0.09) km s -1 arcsec , in agreement with the result of 
1.27 (± 0.10) km s -1 arcsec -1 obtained using the current 
method. 

Regularization can be used to obtain stable solutions 
of Equation (23), but such solutions are not desirable for 
determining the radial behavior of tt p . Regularization 
makes standard tools such as % 2 statistics inappropri- 
ate for calculati ng confidence intervals or evaluating the 
goodness of fit | (|see Aster et al. 20051 ) . This is because 
regularization biases the solution and this bias is often 



a o 




-600 -400 -200 200 400 600 

x (arcsec) 

Fig. 14. — Rotated HI specific intensity map of IC 2574 overlaid 
with contours of constant v v . The galaxy is oriented in the same 
way as shown in Figure 1. The central, most vertical contour is 
v v = km s . Velocity contours are spaced 28 km s" 1 apart. 
The receding half of the galaxy is to the right. The peak flux is 
639 Jy m s _1 beam -1 . In the bottom left corner is a plot of the 
synthesized beam. 

much larger than the confidence intervals that would be 
obtained. Regularization has been tried by Meidt et al. 
(2008b) and M09, but any conclusions based on \ 2 statis- 
tics are unreliable. 

7.2.2. The Solution Method of W98 

When W98 measured the pattern speed of NGC 3031, 
four different, albeit similar, analytic solutions for a con- 
stant flp were derived. Only one such solution is shown 
here to demonstrate the method. It consists of integrat- 
ing Equation (21) once from y = y\ to +00, and then 
again from y = 2/2 > V\ to +00, or similarly for negative 
y. The spacing of y\ and j/2 should be more than a syn- 
thesized beam apart to insure that the calculations are 
independent. After subtracting the two integrations, the 
solution is, 



f_™I"(x,yi)vy(a:,yi)dx - f+™ I„(x,y 2 ) v y (x,y 2 ) dx 

f-™I"fayi)xdx - f+™Iv(x,y 2 )xdx 

(25) 

and similarly for negative y. This solution is equivalent 
to integrating the continuity equation over the region of 
a strip having a width of 1 2/2 1 — \yi\- 

The region of integration in Equation (25) can be much 
smaller than the region in Equation (21). This makes 
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Fig. 15. — Plots of pattern speeds from the solutions for IC 2574. The figure is formatted the same way as Figure 5. Panel (a) shows 
the solution when n = 0. Panel (b) shows the pattern speed having the form of Equation (2) to order n = 1. Panel (c) shows the pattern 
speed having the form of Equation (3) to order n = 1. 
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Fig. 16.— Effect of an incorrect PA for IC 2574. The figure is 
formatted the same way as Figure 6. 



it possible to better constrain the coordinate r in the 
mean value theorem for integration. For each calculation, 
W98 estimated f as the intensity weighted mean radius 
in a strip. Constant and sloping trend lines were then 
fit to the results plotted versus f. The sloping trend 
line provided a statistically significant improvement to 
xt, and a measure of the mean shear rate for the spiral 
pattern of NGC 3031. 

The method of W98 was applied to the H I spiral pat- 
tern of NGC 1365 to see how the results compare with 
those obtained using the current method. It was applied 
over the same region as the current method (\y\ < 400"), 
and used a spacing for y\ and yi of 5 pixels. Individual 
calculations of Equation (25) were spaced 5 pixels apart. 

Figure 18 shows the fits that were made to the results 
using Equations (2) and (3) as forms of f2 p . Due to the 
nature of f, there are no results inside of 170". The un- 
certainties propagated through the calculations of f were 
negligible, therefore the standard errors were found us- 
ing the covariance matrix. There is hardly any noticeable 
difference between panels (b) and (c) of Figure 18. The 
pattern speed in panel (e) appears to be over-fitting the 
results. 

The coefficients for the fitted pattern speeds are shown 
in Table 11. The mean pattern speed and the i = 1 
coefficient for the fit using Equation (3) to order n = 
1 are in agreement with the values obtained using the 
current method. The mean shear rate is exactly the same 
as the value obtained using the current method. 

The statistics for the fits are shown in Table 12. The 
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Fig. 17. — Plots of trend lines fit to the results from applying the original TW84 solution method to NGC 1365. The points and open 
circles represent calculations of Equation (21) for positive and negative y respectively. The solid lines are trends fit to the results. The 
dashed lines are 95% confidence bands. In panel (a) is a constant trend line. In panel (b) is a sloping trend line. 



TABLE 11 

Best-Fit Coefficients for NGC 1365 Using the Method of W98 







ao ± HWo 


ai ± HWi 


a 2 ± HW 2 


Equation 


n 


(km s _1 arcsec -1 ) 


10~ 2 (km s _1 arcsec -2 ) 


10~ 4 (km s _1 arcsec -3 ) 


(2), (3) 





1.07 ± 0.17 






(2) 


1 


2.43 ± 0.53 


-0.55 ± 0.21 




2 


2.25 ± 2.66 


-0.40 ± 2.07 


-0.03 ± 0.38 






10 _1 (km s _1 arcsec -1 ) 


10 2 (km s" 1 ) 


10 3 (km s -1 arcsec) 


(3) 


1 


-4.05 ± 6.00 


3.44 ± 1.36 




2 


-23.1 ± 29.5 


13.1 ± 14.7 


-116 ± 176 



TABLE 12 

Statistics for NGC 1365 Using the 
Method of W98 











P 


R 2 


Equation 


n 


xl 


F 


(%) 


(%) 


(2), (3) 





25.5 








(2) 


1 


13.0 


45.6 


< 0.01 


38.2 


2 


13.3 






38.3 


(3) 


1 


3.5 


289 


< 0.01 


36.1 


2 


0.37 


383 


< 0.01 


38.5 



method of W98 also produces convincing evidence for 
shear in the pattern, and an approximate 1/r behavior 
for fip. The improvement in xl f° r the fit when n = 1 
is extremely statistically significant for both forms of fi p 
(P < 0.01%). For fits using Equation (2), adding the n 
= 2 term did not provide an improvement in xl- For fits 
using Equation (3), adding the n = 2 term produced a 
very small value of xl (= 0.37), which is a sign of over- 
fitting the results. The large values of xl f° r the other 
fits, and the small values of R 2 for all of the radially 
dependent fits, are due to the large amount of scatter 



seen in Figure 18. 

Figure 19 shows plots of the results that include pos- 
sible locations for resonances. When compared to the 
plots in Figure 5, they show how well the results agree 
with those obtained using the current method. The plots 
in panels (a), (b) and (d) resemble those in panels (a), 
(b) and (e) respectively of Figure 5. Note that in panel 
(d) of Figure 19, the fit using Equation (3) to order n 
— 1 is very similar to f2. There are also no clear indi- 
cations of resonances associated with this solution. The 
excellent agreement of the two different sets of results 
confirms the validity of the two different solution meth- 
ods and the robustness of the conclusions arrived at for 
NGC 1365. 

8. COMPARISON WITH PREVIOUS ESTIMATES 
OF n P 

The spiral arm pattern speed of NGC 1365 has been 
estimated previously using other methods that relied 
on the assumptions of the density-wave interpretation. 
JvM95 estimated Q p by identifying possible signs of coro- 
tation and Lindblad resonances in the morphology and 
kinematics of NGC 1365. P. A. Lindblad et al. (1996) 
estimated f2 p by matching barred spiral galaxy models 
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Fig. 18. — Plots of pattern speeds fit to the results from applying the solution method of W98 to NGC 1365. The points represent 
calculations of Equation (25). The error bars represent the uncertainties propagated through the calculations. The solid lines are pattern 
speeds fit to the results. The dashed lines are 95% confidence bands. Panel (a) shows the mean pattern speed. Panels (b) and (c) show 
pattern speeds having the form of Equation (2) to order n = 1 and 2 respectively. Panels (d) and (e) show pattern speeds having the form 
of Equation (3) to order n = 1 and 2 respectively. 



to observations. The pattern speed of a model was set 
by an assumed corotation radius just outside the end of 
the bar. Vera-Villamizar et al. (2001) estimated O p by 
comparing infrared and blue images. They assumed that 
an azimuthal profile phase difference of zero for the two 
images is an indication of corotation. 

The previous estimates of tt p and their corresponding 
radii of corotation resonances are summarized in Table 
13. Originally reported in common units, the previous 
estimates were converted to instrumental units using the 
distances to NGC 1365 that were adopted by the refer- 
ences in the table. Also included in the table are the 
average of the previous estimates and the mean pattern 
speed obtained using the current method. For consis- 
tency with the previous estimates, the uncertainty shown 



for the mean pattern speed is the la uncertainty, which 
is approximated as the standard error. 

The previous estimates are inconsistent with the re- 
sults obtained using the current method. They differ 
from the current result for the mean by at least 3cr, and 
their average is 73% larger than the current result for 
the mean. Furthermore, none of the corotation radii in 
Table 13 are clearly indicated in Figure 5. Also note that 
none of the previous estimates are in agreement with each 
other. They differ among themselves by at least 3cr. 

9. RELEVANCE TO THEORIES OF SPIRAL 
STRUCTURE 

The results for NGC 1365 are inconsistent with funda- 
mental assumptions in the density-wave interpretation 
of spiral structure that are often used for measuring spi- 
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Fig. 19. — Plots of the results for NGC 1365 obtained using the solution method of W98. The figure is formatted the same way as Figure 
5. Panel (a) shows the mean pattern speed. Panels (b) and (c) show pattern speeds having the form of Equation (2) to order n = 1 and 2 
respectively. Panels (d) and (e) show pattern speeds having the form of Equation (3) to order n = 1 and 2 respectively. 



TABLE 13 
Previous Estimates of CI. 





Clp i cr 


Corotation Radius 


Reference 


(km s _1 arcsec -1 ) 


(arcsec) 


JvM95 


2.1 ± 0.1 


138 


P. A. B. Lindblad ct al. (1996) 


1.7 ± 0.1 


157 


Vera-Villamizar ct al. (2001) 


2.7 ± 0.2 


113 


Average 


2.2 ± 0.2 


144 


This work (mean of Cl p ) 


1.27 ± 0.05 
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ral arm pattern speeds, i.e., a constant pattern speed 
with unique locations of corotation and Lindblad reso- 
nances. They are more consistent with the sheared grav- 
itational instabilities of Goldreich & Lynden-Bell (1965), 
the swing amplifier of Toomre (1981), and the recur- 
ring transient patterns in the simulations described by 
Sellwood (2000, 2011). The large amount of shear is 
especially problematic for applying the quasi-stationary 
spiral structure (QSSS) hypothesis of Lin & Shu (1964, 
1966), which proposes that grand-design spiral pat- 
terns are density waves that rotate approximately rigidly 
within a few dynamical timescales. This hypothesis is ap- 
pealing if one assumes that such patterns are long lived 
because it avoids the winding dilemma. For a detailed 
discussion of spiral density waves, its successes in qual- 
itatively explaining many observed properties of spiral 
galaxies, and the QSSS hypothesis, see Bertin & Lin 
(1996) and Bertin (2000, Chapter 18). 

The spiral pattern of NGC 1365 is winding on a char- 
acteristic timescale of ~ 500 Myrs (using the convention 
in M06). A similar timescale was found by M06 for NGC 
1068, and follows from the results of W98 for NGC 3031. 
Winding is regularly seen in simulations of spiral galax- 
ies. Examples can be found in the simulations of isolated 
galaxies by Berman et al. (1978), Bottema (2003), and 
Dobbs & Bonnell (2008); and of patterns generated from 
a tidal perturbation by Sorensen (1985), S. Oh (2008), 
Dobbs et al. (2010), and Struck et al. (2011, in press). It 
is worth noting that the latter two authors provide plots 
of fl p for their simulated galaxies that show an approxi- 
mately 1/r behavior. 

For such a short winding timescale the winding 
dilemma is avoidable if the spiral pattern is a recurring 
transient feature. Support for the existence of transient 
spiral patterns can be found in the redistribution of an- 
gular momentum by spiral arm torqueing (Gnedin et al. 
1995; Binney & Tremaine 2008, Section 6.1.5; Foyle et 
al. 2010), the damping of spiral arm amplitudes by vis- 
cous gas (Toomre 1969), and the inability of simulations 
to produce long-lived spiral patterns (Sellwood 2011). In 
the case of NGC 1365, the bar could regenerate the spiral 
pattern (Sanders & Huntley 1976; Huntley et al. 1978). 

A successful theory of spiral structure must be able to 
account for the shear observed in the spiral arms of NGC 
1365. A more detailed discussion of the consequences 
of the results for theories of spiral structure would be 
premature due to having only applied a rigorous analysis 
to the results for NGC 1365. The new solution method 
should be applied to a sample of galaxies and the results 
thoroughly analyzed before it would be appropriate to 
do so. A preliminary analysis of the results for NGC 
2403, NGC 2903, NGC 3031, NGC 3627, NGC 4535, 
and NGC 5194, however, do show convincing evidence for 
shear in their spiral patterns as well. These results, their 
thorough analysis, and additional tests of the method 
will be presented in forthcoming papers. 

10. SUMMARY 

In this paper a new method was developed for solv- 
ing the TW84 equations. It uses least squares to solve 



an integral equation for simple functional forms of Sl p . 
Standard statistical tools can then be used to determine 
the functional form of fi p that provides the best solution, 
and thus the best estimate of Q p . The method was ap- 
plied to the H I spiral pattern of the barred grand-design 
galaxy NGC 1365. The findings are as follows: 

1. ) The mean pattern speed is 14.5 (± 2.16) km s -1 
kpc -1 . 

2. ) The mean shear rate is - 0.71 (± 0.24) km s -1 
kpc -2 . 

3. ) The best estimate of tt p is 2.87 (± 0.61) xlO 2 km 
s- 1 It. 

4. ) The radial behavior of Q p is approximately 1/r. 

5. ) There arc no clear indications of unique corotation 
or Lindblad resonances. 

6. ) Tests of the method showed that these findings are 
not selection biased. Consistent results were obtained for 
different regions of the galaxy. 

7. ) The pattern speed is very similar to ft. If the 
pattern is decoupled from the material, the difference 
may be too small to distinguish with the precision of 
the results. It is also possible that the assumed PA is 
incorrect, in which case the pattern is most likely rotating 
more slowly than the material. 

8. ) Other methods of solving the TW84 equations for 
shearing patterns were found to produce results in agree- 
ment with those obtained using the current method. 

9. ) Previous estimates that relied on the assumptions 
of the density-wave interpretation are inconsistent with 
the results obtained using the current method. 

10. ) The results are inconsistent with fundamental as- 
sumptions in the density- wave interpretation that are of- 
ten used for finding spiral arm pattern speeds. The re- 
sults are more consistent with spiral structure theories 
that allow for shearing patterns. 

11. ) The characteristic timescale for winding is ~ 500 
Myrs. The winding dilemma is avoidable if the spiral 
pattern is a recurring transient feature. 
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